Analytical formulas for gravitational lensing: higher order calculation 
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I. INTRODUCTION 

In a recent paper, [l] we have introduced a new method for the calculation of the deflection angle of light in a 
general static and spherically symmetric metric in general relativity. This approach allows one to convert the integral 
for the angle into a geometrically convergent series, whose terms can be analytically calculated. Because of the strong 
rate of convergence of the series, typically few terms are already sufficient to obtain highly accurate analytical results. 
For this reason, previous analysis |l| was focussed on the calculation of first order, and a general formula, valid for 
an arbitrary static and spherically symmetric metric, was derived. The series is not a typical perturbative expansion 
in some small physical parameter, and therefore is valid also in proximity of the "photon sphere" . 

A good description of the deflection angle close to the photon sphere is particularly important for the study of strong 
gravitational lensing, a topic which has received wide attention in the recent past: for example strong gravitational 
lensing in a Schwarzschid black hole has been considered by Frittelli, Kling and Newman Q and by Virbhadra and 
Ellis Virbhadra and collaborators have also treated the strong gravitational lensing by naked singularities 
and in the presence of a scalar field [j|; Eiroa, Romero and Torres |g| have described Reissner-Nordstrom black hole 
lensing, while Bhadra has considered the gravitational lensing due to the GMGHS charged black hole 0; Bozza 
has studied the gravitational lensing by a spinning black hole jg ; Whisker [|| and Eiroa |lfjj have considered strong 
gravitational lensing by a braneworld black hole; still Eiroa ^lj has recently considered the gravitational lensing by an 
Einstcin-Born-Infcld black hole; Sarkar and Bhadra have studied the strong gravitational lensing in the Brans-Dickc 
theory [T^: finally Perlick |l3j has obtained an exact gravitational lens equation in a spherically symmetric and static 
spacetime and used to study lensing by a Barriola-Vilcnkin monopolc and by an Ellis wormhole. 

To describe the strong gravitational lensing regime Bozza has introduced and analytical method based on a careful 
description of the logarithmic divergence of the deflection angle which allows to discriminate among different types 
of black holes, but which can only be used in proximity of the photon sphere^ij]. O n the other hand, other methods, 
which are very precise in the weak lensing regime have also been developed: for example, Mutka and Mahonen [l5lll6|| 
and Belorobodov 01 have derived improved formulas for the deflection angle in a Schwarzschild metric, while Keeton 
and Petters0 have introduced a formalism for computing corrections to lensing obscrvablcs in a static and spherically 
symmetric metric beyond the weak deflection limit. These methods fail to describe the strong lensing regime and do 
not even reproduce a photon sphere. This is in contrast with our method of Q, which already to first order predicts 
a photon sphere and provides much better results than usual perturbative expansions. 

The purpose of this paper is to further improve the description of the strong lensing regime extending previous 
analysis [I| to higher orders and providing with a general formula for the systematic and analytical calculation of the 
deflection angle in any static spherically symmetric metric. 

The paper is organized as follows: in scction[n]we review our method. In scction lTTTl we obtain the general expression 
for the higher order contributions. In section ITVl we discuss the application of the series to a few examples. Finally in 
section IVI we draw our conclusions. 
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II. THE METHOD 

In this section we outline the method introduced earlier Q| .We consider a general static and spherically symmetric 
metric which corresponds to the line element 

ds 2 = B{r)dt 2 - A{r)dr 2 - D{r)r 2 (d6 2 + sin 2 9 d(f> 2 ) . (1) 

Notice that the Schwarzschild metric is a special case of this equation. Although at this stage we leave the metric 
coefficients A(r), B{r) and D(r) unspecified, we assume that linv^oo f(r) = 1, where f(r) = (A(r), B(r), D(r)) which 
are consistent with a flat spacetime at infinity. 

One can express the angle of deflection of light propagating in this metric in terms of the integral [l^ 



A0 = 2 / y/A(r)/D(r) 



2 D(r) B(r ) "* 
r J D(r ) B(r) 



( 2 ) 



where ro is the distance of closest approach of the light to the center of the gravitational attraction. 

In general it is possible to evaluate analytically Eq. only in a limited number of cases, depending upon the form 
of the metric. Nonetheless we wish to prove that precise analytical formulas can still be obtained even though no 
explicit expression for the exact integral may exist. 

As explained in Q] one can introduce the function 

2 D(r Q /z) D 2 (r /z) B(r Q ) B(r ) 

[ 1 ~ A(r /z) A(r /z)B(r /z)D(ro) + D(r ) ' 1 1 

where z = r$/r and thus convert Eq. J5J) into the form 

A0 = 2 /' - dz -tt . (4) 

Jo y/V(l)-V( Z ) 

This expression bears a close resemblance to the expression for the period of a classical oscillator with a potential 
V(z). It was shown in Q that an efficient way of dealing with such integral is to use the nonperturbative method 
based on the Linear Delta Expansion (LDE), developed by Amore and collaborators [20ll2l| . 

Let us review the main aspects of the method: we interpolate the full potential V(z) with a simpler potential 
Vq(z), which should be chosen in such a way that the integral inside Eq. |@} can be performed explicitly when V(z) 
is substituted by Vq(z). Then, we write 

V s (z) = V (z) + S(V(z)-Vo(z)) , 

where 6 is a dummy expansion parameter that is chosen equal to unity at the end of the calculation. The reference 
potential Vq(z) may depend upon a set of adjustable parameters, which we will collectively call A. It is sufficient for 
present purposes to consider the simplest reference potential Vq(z) = Xz 2 and rewrite the deflection angle as 

„ f 1 dz I 
A0 5 = 2/ — ===— === -7r, (5) 



where E — Vb(l) = A, 



o ^/Ev-Vo(z) y/l + SA(z) 



AW,fa-l. (6) 

Eq - Vo(z) 



and E = V(l). Clearly Eq. JSJ reduces to the exact expression Eq. 10} when 5 = 1. 

If |A(z)| < I for < z < 1 the expansion of equation JSJ) in p owers of 5 yields a perturbation series that converges 
towards the exact result. It was shown in earlier papers [l|,|20j|21| that this condition is met when A is greater than 
a critical value, A > Ac- Thus, the sequence of partial sums converges towards the exact result that is independent 
A, although each partial sum depends on the adjustable parameter. In order to minimize the dependence of the 
approximate results on A we resort to the the Principle of Minimal Sensitivity (PMS)[22| which in present case reads: 



(7) 
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where Acp^ is the partial sum of order N. 

This condition determines an optimal value of A around which the truncated series is less sensitive to changes in 
that parameter. In addition to it, one expects that the resulting series exhibits improved convergence properties @|. 
More precisely, it has been verified that the series obtained by means of the PMS exhibit geometric rate of convergence 
providing markedly accurate results with just a few terms In fact, just a first order calculation proved sufficient 
for the prediction of the location of the photon sphere with remarkable accuracy Q . 

Previous calculations pj also show that the rate of convergence of the PMS series depends critically upon the 
physical parameters of the model, and in particular the accuracy of the results deteriorates as tq approaches the 
photon sphere. For this reason in this paper we are interested in the systematic improvement of our results by means 
of general and simple formulas that facilitate the calculation of any desired order of approximation for a given arbitrary 
static and spherically symmetric metric. 



III. HIGHER ORDER CALCULATION 



Under the assumption that the general metric considered here is flat at the infinity we can express the potential in 
Eq. J2Jl as 



oo 



V( Z ) = J2vnZ n . (8) 

Substitution of this expression into Eq. J5J) yields 

oo n— 1 



71 = 1 



n=l k=Q 

Additionally, the angle of deflection becomes 

A0 = 2 C . dz y (2n ~f )!! S»(-l) n A n {z) -,eVA„, (10) 

where the meaning of the terms A„ is obvious. 

Assuming that |A(z)| < 1 for z € (0,1) and the series is convergent, then summation and integration can be 
interchanged, thus leading to 

^^EE 2 i^.S S " £ V ri V r2 ...V rj 0« 2 ... r .-7T (11) 
ri=0j=0 ^ ' r± r 2 ... rj 

where the rank- A tensors 

n (N) \ \ V / - (12) 

"nin 2 ...njv — L^l L (l + z) N Vl - r2 ' 



fc 1= k 2 =0 k N =0 



are dimensionless and completely symmetric in the lower indices. The explicit form of these tensors up to order four 
is given in the Appendix. 

By induction one derives the general expression 



r>A:-l 



^_ iXE|LWg+Vg) y 

-l)!! V?r | r(E- =1 n,/2-(fc-l)) 



r(EU-"*/2 + i/2) 



(2k 1)!! 1 r(Ef =1 nj2 - (k - 1)) ^ r(E-^ nj2 (k - 1)) 



+ E F( ^" /2 + 1/2) + - ■ (13) 
where p means sum of all the different combinations of the k — I elements chosen out of k total elements. 

31 32 ■ ■ H 

The amount of terms is simply given by the binomial coefficient ( , ^ , j = u(£!_iy m ■ 
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A crucial observation is that the indices of f2 take all possible values when the sum is performed. Consequently 
we cannot distinguish among the possible different combinations inside 113fl and the same result must be obtained 
taking a single one with the proper multiplicity factor given by the binomial. We can therefore make the substitution 
o(fc) spk A(fc|0 where 

nl - n,_ (2k-i)\\ Vw {k-i) ( 1} r(ESi/2-(fc-i)) " ( ' 

By substitution of this expression into Eq. i|llf> we obtain the final expression 

^tt ^ ( ~3^ Sn t ^v r2 ... Vrj ±nBn-^ (15) 

n=0j'=0 J ' V J '' n r 2 ... fj 1=1 

In the end, we can explicitly write the expression corresponding to the first few orders as 
\a m 315 y 189 45 35 

A0 = -7T + , S ^ Si H ry^— S 2 S 3 H jyj— £ 4 + ... (16) 

64VA^ 16A 3 P ^ S 32A 5 P ^ S 16A^ S UX P /2 MS 

where 

oo j 

Ej= J2 V T1 V r2 ...V rj ^2(i j S, ■ (17) 
n r 2 ... rj ;=i 

For simplicity, in Eq. (|16[) we have selected the optimal parameter of first order A = Xpms for all orders of 
approximation. The reason is that it is not possible to obtain analytical expressions of Xpms for greater orders, and 
that it is sufficient that Xpms > Xc for the series to be geometrically convergent. Thus, in all the applications below 
we choose the value of A pms obtained previously 0] : 

(1) gi_ 2 v r(n/2 + l/2) 
XpMS ~ V^ 1 - T(n/2) ■ (18) 



IV. APPLICATIONS 



In this section we consider some applications of the general formula l|15|) . The first example is the Schwarzschild 
metric given by 



B(r) = A~ 1 {r) = (l- 



2GM\ 



, D(r) = 1 , 



(19) 



where M is the Schwarzschild mass and G the gravitational constant. The angle of deflection of a ray of light reaching 
a minimal distance ro from the black hole is given Eq. J3J , and the exact result can be expressed in terms of incomplete 
elliptic integrals of the first kindj^ as 



A<£ = 4W-^ F 1—,k) — F ((p, k) 



(20) 



where tq = rO/GM and 



T EE 



r -2 
r + 6 



ee y(T^r+6)72T , 



arcsm 



2 + T - r o 



6 + T - r o 



(21) 



The "potential" for this metric is particularly simple since it is just a polynomial of third degree: 



V{z) = z 2 



2GM 

ro 



(22) 
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FIG. 1: Error of the deflection angle as a function of ro for the Schwartzchild metric with GM = 1 . The dotted vertical line 
shows the location of the photon sphere, (color online) 



In order to apply Eq. (|15|l we choose GM = 1 and calculate the error of the deflection angle as a function of the 
distance of maximal approach ro (= ro), defined as 



A(j)pMS - A(/> e 



A0 e 



x 100 , 



(23) 



where A<fi exact is the Darwin solution (|2U|) . Fig. shows that the errors of the PMS approximations through fourth 
order arc extremely small and smaller that one percent even close to the photon sphere located at ro = 3 (the one 
percent error is obtained at ro ~ 3.2 for the fourth-order partial sum). 

Present approach yields similar results for the Reissner-Nordstrom (RN) metric 



2GM Q 2 

B{r)=A-\r)= [ 1 - — — + ^ 



^ > D ( r ) = 1 



that describes a black hole with charge. In this case Eiroa, Romero and Torres | 
deflection angle in terms of elliptic integrals of the first kind (see eqn. (A3) of 0). 
The "potential" for this example is also polynomial: 



V{z) 



2GM 



Q 2 



(24) 

have been able to express the 



(25) 



The behavior of the error displayed in Fig. [21 is similar to that observed in the case of the Schwarzschild metric. 

Both examples discussed above lead to polynomial potentials in z and are amenable to exact solutions in terms of 
elliptic functions. However, it should be noticed that our method applies to more general cases that cannot be solved 
exactly. 

In an earlier paper 0] we considered the propagation of light in a charged black hole coupled to Born-Infeld 
electrodynamics |ll). which corresponds to the effective metric 



A(r) 



i>(r) 



B(r) = vMOV-M , D{r) = 



where 




F 



\bQ\ 



■ 2 + \bQ\) ' V2 



(26) 

(27) 
(28) 
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FIG. 2: Same as in Fig. Q for the RN metric with GM = 1 and Q = 1. (color online) 



TABLE I: Deflection angle for the Born-Infeld metric with GM = 1, Q = 1/2 and 6=1. 



ro Numerical PMSi PMS 2 PMS Z PMS± 

3 5.04066558 3.9836332 4.769690116 4.695977905 4.961794995 

4 1.989492408 1.905459664 1.988326336 1.985415327 1.989578734 

5 1.302767649 1.274857482 1.302841883 1.302239811 1.302781935 
10 0.4908767995 0.4881229478 0.4908905773 0.4908704512 0.4908768924 
50 0.08300041129 0.08293664219 0.08300048071 0.08300040744 0.0830004113 
100 0.04073430348 0.04071929339 0.04073431165 0.04073430326 0.04073430348 



Here, F(a,b) is the incomplete elliptic integral of first kind. We follow the convention of Ref. [Tlj and set G = 1. 
In this case one obtains the potential 



V(z) ee ~ 2( 7 V +Tp( r oMro) (29) 

uj{r /z) uj z {r /z) 



which can be expanded around z = as 



V{z) = Y,v n z n . (30) 



Clearly, in this case the accuracy of our results will depend not only upon the order of the partial sum but also 
on the number of Taylor coefficients used to approximate the potential. We illustrate this point in Fig. [3] for the 
deflection angle when GM = 1, Q = 1/2 and 6 = 1, and for ro = 3, which is very close to the photon sphere where 
the angle diverges. This situation clearly corresponds to a strong gravitational lensing regime. The horizontal axis 
shows the number of Taylor coefficients of the potential used in the calculation; the reader will appreciate that just a 
small number of them (n « 4) is sufficient for reasonable accuracy. For larger values of n the results do not appear to 
improve noticeably and a sort of plateau is reached. We point out that the main advantage of our approach is that 
it leads to relatively simple analytical functions of the parameters of the metric. 

Table |U shows the deflection angles in terms of ro for the Born-Infeld metric with the same parameters as before. 
The potential V(z) is approximated with a Taylor polynomial of order 10 and the calculations are performed with 
our method through order 4. 

As a final example we consider the metric of Weyl gravity [24L I2M l26l |2^ | : 



--A- 1 (r)= (l-^-+ 7 r-fcr 2 ) 
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FIG. 3: Deflection angle for the Born Infeld metric with GM = 1, Q = 1/2, 6=1, and ro = 3. The horizontal axis specifies 
the order of the Taylor approximation to the potential. The dashed line is the numerical result, (color online) 



where /3,7 and k are constants. When (3 = GM and 7 and k small enough one recovers the Schwarzschild metric on a 
certain distance scale. The linear and quadratic terms are significant only on galactic and cosmological scales, if the 
corresponding parameters are sufficiently small. 

The potential corresponding to this metric is once again polynomial in z: 

Viz) = -kr 2 + irQ z + z 2 -—z 3 . (32) 

Notice that the constant term cancels out of the integral for the deflection angle that will therefore be independent 
of k. 

On defining A = ^Jl — + 2212. we obtain the fourth-order expression 

A^-r + E^, (33) 



\2n+l ' 

n=0 



where 



di = 


3.19285 






d 3 = 


0.42847/3 


- 0.10350 






ro 






d 5 = 


1.02342/3 2 


0.48610/3 


- 0.058569 




r 2 
'0 


ro 




d 7 = 


0.75158/3 3 


0.55570^ 2 


! 0.13858/3 




r 3 


r 2 
'0 




dg = 


1.15417/3 4 


1.18963^ 3 


0.46370/3 
+ T 2 - 




r 4 
'0 


r 3 





(34) 
(35) 

(36) 

0.01166 (37) 
f 0.00533 . (38) 



ro 

Expanding this equation in powers of (3 and 7 we obtain 

A^)- 4 / 3 ^f^P 2 4/3 2 3^(3 j 2 r 2 \ 

V)«-- 7 n,+ ^- 5 a--- ; 3-+7/3- — + — J+... , (39) 

where the term between parenthesis contains the contribution of second order in the expansion. Notice that the first 
two terms are those usually reported in the literature [25I [2tJ . 
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Fig. 2] shows results for the deflection angle given by the Weyl metric with j3 = 7 = 1. Our fourth-order analytical 
formula agrees remarkably well with the numerical calculation. 

Notice that the deflection angle displays an interesting behavior depending upon the sign of 7: as we see in Fig. 0] 
for 7 > the angle is negative at large distances while it becomes positive close to the photon sphere. Using our 
analytical formula we are able to provide an accurate expression for the value of tq corresponding to absence of 
deflection: 

[p /3(-/3(40960 + 97r(112 + 757r(-16 + 37r)))7 + 864^tt(- 16 + 5ir)^y + 1024(32 - 9tt)) 
r *~ Y 7 + 32768 ' *■ ' 

In the case of Fig.0|this formula predicts the location of r* ~ 2.08, which is remarkably close to the exact numerical 
value. 

In the opposite regime, 7 < 0, the deflection angle is real only on a finite region that we obtain from 

apms — t 1 . (41 j 

2 

In fact, when (3 > and 7 > — then Xpms > only for 



^64/37 + tt 2 7T ^64^7 + tt 2 tt . . 

A T" < ^ < A ~A~ ■ ( 42 ) 

47 47 47 47 

Iii Fig. [SI we compare our fourth-order analytical result with the numerical one for the Weyl metric with (3 = 1 and 
7 = — 7r 2 /128. The agreement is remarkable for all the values of rg considered. 



V. CONCLUSIONS 



In this paper we have generalized a method developed earlier [l| for obtaining analytical expressions for the deflection 
angle of light travelling in an arbitrary static and spherically symmetric metric. We have applied a new higher-order 
formula to several examples confirming the accuracy suggested by a rigorous proof and estimate of the rate of 
convergence jl[. 

We believe that our formalism is particularly useful for several reason: first of all, because it stands on a firm 
mathematical ground; second, because it is completely general and can be applied with little or no effort to any such 
metric; finally because the analytic results obtained with our method could provide a valuable tool for the analysis 
of different models. Just to mention one example, in the case of the Weyl metric, we have been able to derive an 
analytical expression for the absence of deflection for a ray of light travelling in a metric with (3 > and 7 > 0. 
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FIG. 5: Deflection angle for the Weyl metric with (3 = 1 and 7 = — 7r 2 /128. 



APPENDIX A: O 



Ci) 



Explicit expressions of the symmetric tensor f2 through order four, which one obtains by straightforward algebra: 



Q( 2 ) 

nin2 



" nin 2 ri3 



ni n2n3?i4 



-T(n/2+ 1/2) 
r(n/2) 

^/ r r (^ti) , r + n 2 + i))' 



-V7T - 



4 

15 " 



r(f-i) r(f-i) r(i(n 1 + na -2)) 



v 2 



r(^-2) 



r(f - 

r(fgg) 
r(f -2) 



' »3 + l ' 

, 2 , 



r(i(m + n 2 + i)) 



r(f -2 



£(| (m+rt 3 + l)) _ r(| (n 2 +n 3 + l)) 

4)) 

rf 



) r(i(m + n 2 -4)) 

r (| (ni + n 2 + n 3 + 1)) 



r (i (nx + n 3 - 4)) r(i(n 2 +n 3 

8 _ / r 



' rai + l ' 
> 2 , 



v 2 ; 



2 , 



(nx + n 2 + n 3 

p | , n 4 + l 'j 



4)) 



105 v \ 

r(i(m 



ILL 

2 



--3) r( 
1)) , r(i(m 



n 3 



r 
D) 



-3) 
r(i(n 2 



n 3 



-3) 

D) 



r(i(m+n 2 -6)) r(i(m+n 3 -6)) r(±(n 2 + 

r 



(Al) 
(A2) 

(A3) 



(A4) 



"3 - 6)) 

r(i(m + n 4 + i)) r(l(n 2 + n 4 + l)) r(i(». 3 + n 4 + i)) 
r (|(m + n 4 - 6)) + r (i(na + n 4 - 6)) + T (±(n 3 + n 4 - 6)) 
r(|K +n 2 + n 3 + l)) r (l(nx + n 2 + n 4 + 1)) 
T (|(ni + n 2 + n 3 - 6)) T (±(nx + n 2 + n 4 - 6)) 

r(i(nx+n 3 +n 4 + l)) r (±(n 2 + ?i 3 + n 4 + 1)) r (±(ni + ?i 2 + n 3 + n 4 + 1)) ' 
r(i(m +n 3 + n 4 -6)) r ( ±(n 2 + n 3 + n 4 - 6)) r (|(m + n 2 + n 3 + n 4 - 6)) 



(A5) 
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APPENDIX B: E fc 



Explicit expressions for the coefficients defined in Eq. I|17|) for the Taylor expansion of an arbitrary potential 
through order three: 

50 = \ (Bl) 

51 = vi + -TTV 2 + 2v 3 + ... (B2) 

2vj 2 1 2 

5 2 = — + 2V2VI + TTV3VI - -V3 V\ + -TTV 2 

+ ^--3+4^3 + ... (B3) 

5 3 = + 2v 2 v\ + ^-v 3 vj + 3vj vi - \irvlvi 

15 5 2 

47 2 „ 1 3 3 , 

+ ~rv 3 vi + 3ttv 2 v 3 Vi - 2v 2 v 3 vi + -ttv 2 - -n v 3 
5 2 2 

+ + —^v 2 v% - Av 2 vi + 6vp 3 + . . . (B4) 

lSV-i 28 o 76 o - 9 9 9 9 

5 4 = -^- + ^v 2 vl + —v 3 vl + Av 2 2 vl + 3nv 2 3 vl 

124 2 2 48 2 , 3 n 3 

- "35" w 3 w i + y «2«3«i + 4 *>2 w i + 9nv 3 vi 

708 , „ 2 188 2 

- -^- w 3 w i - 67rw 2 w 3 z;i + — v 2 v 3 Vi 

2 2 1 4 99 4 104vf 

+ 6tTW 2 W3U1 - 4v 2 V 3 Wl + -TTV 2 + — 7T 1> 3 — 



488 15 

6irv 2 vl + — —v 2 + —ttv 2 v 2 - 8«2«3 + 8v 2 v 3 + . .. (B5) 

_L A 
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